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ABSTRACT 

On  the  basis  of  some  approximate  equations  derived  for  nonlinear  water 
waves  in  a  channel  of  variable  cross  section,  several  problems  of  interest, 
which  have  attracted  such  attention  recently,  will  be  taken  up  in  this 
report.  First  we  shall  consider  the  breaking  of  an  acceleration  wave  moving 
toward  a  shoreline  in  a  general  channel.  Next  the  so-called  infinite  mass 
dilemma,  which  arises  in  the  study  of  the  development  of  a  solitary  wave  in  a 
channel  of  variable  cross  section,  will  be  resolved.  Finally  we  shall  use  an 
approximation  method  to  study  the  fission  of  solitons  in  a  general  channel  and 
justify  it.  i 
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SIGNIFICANCE  AND  EXPLANATION 


Consider  a  wave  in  an  open  channel  moving  toward  a  shoreline*  Under 
physically  reasonable  assumptions*  we  shall  determine  whether  the  wave  will 
break*  and  the  location  of  the  break  if  it  does.  In  recent  years*  there  has 
also  been  growing  interest  in  the  development  of  a  solitary  wave  in  an  open 
channel  of  variable  cross  section.  It  is  found  that  a  shelf  is  generally 
formed  behind  the  solitary  wave.  In  some  studies*  the  shelf  is  found  to 
extend  to  infinity.  This  is  the  so-called  infinite  mass  dilemma.  In  other 
words*  infinite  mass  would  be  created  or  annulled  by  a  perturbation  on  the 
solitary  wave.  We  shall  resolve  this  dilemma  by  showing  that  the  shelf  can 
only  be  finite.  On  the  other  hand*  if  a  solitary  wave  propagates  from  one 
uniform  cross  section  of  a  channel  to  another*  the  solitary  wave  may  split 
into  several  soli tons.  We  shall  use  an  approximation  method  to  find  a 
criterion  for  the  fission  of  n  solitons  in  a  general  channel  and  show  that 
the  approximation  method  indeed  yields  an  approximate  solution  to  the  exact 
solution  of  the  governing  equation. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


SOME  NONLINEAR  PROBLEMS  OP  MATER  NAVES  IN  A  CHANNEL 


M*  C.  Shen 

1 .  INTRODUCTION 

In  this  report  ms  shall  study  several  interesting  problems  of  nonlinear  water  waves  in 
a  channel  of  variable  cross  section  by  means  of  some  approximate  equations  derived  from  the 
hydrodynaaical  equations*  For  finite  amplitude  waves,  a  system  of  shallow  water  equations 
will  be  given*  For  small  amplitude  waves,  we  shall  use  a  K-dV  equation  with  variable 
coefficients  derived  in  a  previous  report  (Zhong  and  Shen,  1982). 

One  of  the  fascinating  phenomena  of  finite  amplitude  water  waves  in  a  channel  concerns 
the  breaking  of  a  wave  moving  toward  a  shoreline,  the  development  of  a  bore  and  the 
movement  of  the  shoreline  after  the  bore  reaches  it.  For  the  two-dimensional  case 
corresponding  to  a  rectangular  channel  of  variable  depth,  the  bore  run-up  problem  was 
studied  by  Roller,  Levine  and  Whithaa  (1960),  Ho  and  Meyer  (1962),  and  Shen  and  Meyer 
( 1963a, b)  on  the  basis  of  shallow  water  equations  (Stoker,  1957).  Later  Gurtin  (1975) 
derived  a  criterion  for  the  breaking  of  an  acceleration  wave  in  a  two-dimensional  channel, 
and  his  result  was  extended  by  Jeffrey  and  Mvungi  (1980)  to  the  case  of  a  rectangular 
channel  of  variable  width  and  depth  Me  shall  make  use  of  the  shallow  water  equations  for 
a  general  channel  and  generalise  Gurtin* s  result  to  predict  the  breaking  point  of  an 
acceleration  wave*  Needless  to  say,  the  use  of  shallow-water  equations  for  the  study  of 
bore  propagation  may  be  open  to  criticism*  The  issue  would  be  settled  if  we  knew  the 
precise  conditions  for  the  validity  of  these  approximate  equations*  Up  to  date,  the 
shallow  water  equations  for  a  two-dimensional  channel  with  analytical  initial  data  have 
been  justified  by  Kano  and  Nishida  (1979),  and  for  a  nearly  uniform  channel  with  a  priori 
assusq>tlons  on  the  free  surface  by  Berger  (1976).  At  present  we  may  accept  shallow  water 
equations  as  model  equations, and  the  bore  run-up  problem  for  a  general  channel  certainly 
deserves  further  investigation. 


Sponsored  by  the  United  States  Army  under  Contract  No*  DAAG29-80-C-0041  and  the  National 
Science  Foundation  under  Grant  No.  MCS800-1960. 


Recently  the  problem  of  the  development  of  «  solitary  wave  in  a  channel  of  variable 
cross  section  has  attracted  much  attention,  and  there  have  been  discussions  on  the  so- 
called  Infinite  mass  dilemma,  which  arises  from  the  formation  of  a  shelf  behind  the 
solitary  wave.  If  the  shelf  were  extended  to  infinity,  then  infinite  mass  would  be  created 
or  annulled  by  a  perturbation  on  the  solitary  wave.  A  study  of  this  problem  may  be  found 
in  Miles  (1979)  and  Knickbocker  and  Newell  (1980),  based  upon  the  K-dV  equation  for  a 
rectangular  channel  of  variable  depth  or  width  (Kakutanl,  1971;  Johnson,  1973;  Shuto, 

1974).  As  an  extension  of  the  existence  results  due  to  Kato  (1975,  1980),  a  global 
existence  theorem  will  be  established  here  for  the  solution  of  the  K-dV  equation  with 
variable  coefficients.  It  follows  that  the  shelf  can  only  be  finite  if  it  is  formed  behind 
the  solitary  wave. 

To  study  the  fine  structure  of  the  shelf  behind  a  solitary  wave,  we  shall  consider  the 
solitary  wave  moving  from  one  uniform  cross  section  of  a  channel  to  another  through  a 
transition  section.  This  problem,  in  fact,  has  a  long  history.  Using  a  system  of 
Boussinesq  equations,  Madsen  and  Mel  (1969)  studied  numerically  the  problem  of  a  solitary 
wave  propagating  from  a  channel  of  constant  depth  to  a  shoal  of  constant  but  smaller  depth 
through  a  transition  region  with  mild  slope.  They  found  that  the  solitary  wave  is 
disintegrated  into  a  train  of  solitons  of  decreasing  amplitudes,  In  qualitative  agreement 
with  experiments.  An  analytical  investigation  of  the  shoal-induced  fission  of  solitons  was 
made  by  Tappert  and  Zabusky  (1971).  They  used  the  WKB  approximation  to  describe  the 
transition  of  the  solitary  wave  and  the  inverse  scattering  method  for  the  K-dV  equation 
with  constant  coefficients  by  Muira  et  al.  (1974),  and  obtained  a  criterion  for  the  fission 
of  n  solitons.  Johnson  (1973)  derived  a  K-dV  equation  with  variable  coefficients  for  a 
channel  of  variable  depth  and  found  the  same  criterion  by  a  similar  method.  Recently 
Djordjevic  and  Redekopp  (1978)  considered  the  disintegration  of  internal  solitary  waves  and 
Zhong  and  Shen  (1982)  extended  the  previous  results  to  the  case  of  a  symmetric  triangular 
channel.  As  a  simple  consequence  of  the  existence  results  just  obtained,  the  approximation 
method  used  is  shown  to  be  valid  if  the  distance  to  cover  the  transition  section  is 


sufficiently  small 


In  Section  2,  we  derive  the  shallow  eater  equations ,  and  study  the  breaking  of  a  wave 
in  a  general  channel*  In  Section  3  we  establish  the  global  existence  theorem  for  the  K-dV 
equation  with  variable  coefficients  and  resolve  the  infinite  mass  di leave*  Finally  the 
approximation  method  is  used  and  justified  for  the  study  of  fission  of  solitons  in  a 
general  channel  in  Section  4* 


2*  SHALLOW  MATER  SQUAT  IONS  AND  THE  BREAKING  OF  A  WAVE 


/***a*KsaM 


Nt  consider  the  ir rot at ion* I  Motion  of  an  inviscid,  incompressible  fluid  of  constant 
density  under  gravity  in  a  channel  with  a  boundary  defined  by  h*(x*,y##**>  ■  0/  where 
s*  is  positive  upward  and  x*  is  in  the  longitudinal  direction  (Fig.  1).  The  governing 
equations  are 

V*-q*  -  0  ,  (2.1) 

V*xq*  -  0  ,  (2.2) 

p(q£#  +  q# • Vq# )  -  -Vp*  +  g  ,  (2.3) 

subject  to  the  boundary  conditions 

♦  q*»V*z*  -  0,  at  z*  -  -a*  +  n#(t*,x*,y*)  -  0  ,  (2.4) 

p*  -  0  ,  (2.5) 

q#*V*h#  -  0,  at  h*  *  0  •  (2.6) 

Here  V*  ■  ( 3/3x*, 3/3y#, 3/3 z#),  q#  -  (u*,v#,w*>  is  the  velocity,  t#  is  the  time, 
g  *  (0,0, -g)  is  the  constant  gravitational  acceleration,  p  is  the  constant  density, 
p*  is  the  pressure,  and  s#  *  n*  is  the  equation  of  the  free  surface.  To  derive  the 
shallow  water  equations,  we  make  the  following  assumptions.  The  channel  boundary  is 
convex,  sufficiently  smooth,  and  varies  slowly  in  the  longitudinal  direction;  the  magnitude 
of  the  transverse  velocities  is  much  smaller  than  that  of  the  longitudinal  velocity.  In  a 
manner  as  suggested  by  Friedrichs  (1948),  we  introduce  nondiraenslonal  variables 
t  -  0-V2tV(H/g)1/2,  (x,y,r)  -  ( 0_1/2x*/H,y*/H,*VH)  , 

n  -  n*/H,  h  -  h*/H,  (u,v,w)  -  (uV(gH)1/2.81/2v*/(gH)1/2,01/2w/(gH)1/2)  , 

1/2 

where  0  *  L/H  and  L  and  H  are  respectively  the  horizontal  and  transverse  length 

scales.  In  terms  of  them,  (2.1)  to  (2.6)  become 


u  +  v  ♦  w  •  0  , 
x  y  z 


0u  -  V  • 
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v.  ♦  uv  + 
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wv  +  0p  -  0  , 

t  X 

y 

z  *y 

+  uw  +  vw 
x  y 

♦  ww 

z 

♦  0(pz  +  1)  -  0 

(2.7) 

(2.8) 

(2.9) 

(2.10) 

(2.11) 


n  ♦  un  ♦  vn  -  w  -  o 
t  x  y 


at  z  ■  n 


(2.12) 


(2.13) 


uh  +  vh  ♦  xh  *0  at  h  ■  0  *  (2.14) 

x  y  z 

Aim  that  u,  v,  w,  p  and  n  possess  an  asy^totic  expansion  of  the  for* 


*  ~  *0  +  6-1  *1  +  0-2 *2  +  *•*  ' 


(2.15) 


and  substitute  (2.15)  in  (2.7)  to  (2# 14).  The  equations  for  the  zeroth  approximation  are 


uA  ♦  vA  +  *  0  . 

Ox  Oy  Os 

uA  *  *  0  , 

Oy  Oz 

U0t  +  U0u0x  +  Pox  +  V0u0y  +  W0u0z  "  0 

Poy  “  °»  Pox  “  _t  ' 

”ot  ♦  Vox  +  voV  -  *0  *  0  1  . 


„  -  s.  •  0  ) 

)y  0  I 

Pq  *  0  / 


(2.16) 

<2.17) 

(2.18) 

(2.19) 


(2.20) 


(2.21) 

(2.22) 


Ughjj  +  vQhy  +  e0hE  *  0  at  h  ■  0  .  (2.22) 

As  seen  frosi  (2.17),  (2.19)  and  (2.21 ),  Uq  is  a  function  of  t,x  only,  and 

P0  *  -s  ♦  ,  (2.23) 

which  implies  nn  is  also  a  function  of  t#x  only.  It  follows  from  (2.17),  (2.18)  and 


(2.23)  that 


u_^  4-  u  u„  +  rv  *  0  . 
Ot  0  Ox  T)x 


(2.24) 


Now  we  integrate  (2.16)  over  a  cross  section  D  of  the  channel,  apply  the  divergence 
theorem  and  make  use  of  (2.20)  and  (2.22)  to  obtain 


//  (v0y  +  *0z,dyd*  “  -tt0xMt'X)  * 
D  * 


•  -  “o  /  hx<hy  +  *  (n0t  +  UoV,B<t'X>  ' 

I» 


By  rearranging  the  terms,  we  have 


n0t  +  U0n0x  +  u0xA(t,x)/B(t,x)  -  (u0/8(t,x)l  /  h^(h2  ♦  h2)~1/f2da  -  0  , 


(2.25) 


x  g  * 


where  A(t,x)  is  the  area,  B(t,x)  is  the  width  and  L  is  the  wetted  boundary,  of  the 
cross  section  D  (Fig.  1).  (2.24)  and  (2.25)  form  a  system  of  nonlinear  equations,  which 


7  * 


V  .*  A  -A  -A  *  .  A  .’ 


uy  bt  uMd  to  model  boro  foraatlon  «nd  its  subsequent  dmlogmt  in  a  channel  of  variable 
cross  section. 

In  the  following  we  extend  Curtin**  Method  to  the  case  of  a  general  channel*  The 


assumptions  made  are  the  following i 

( 1)  uQ,  are  continuous# 

(2)  the  first  and  second  derivatives  of  Uq  and  nQ  possess  at  most  jump 
di scontinui ties # 

(3)  «0  *  ■  0  ahead  of  the  wave* 

Denote  the  value  of  a  function  f  immediately  behind  the  wave  front  by  f~*  Hereafter  we 


also  drop  the  subscript  0*  From  assumptions  (1)#  (2)v  we  have 

u  ■  n  ■  o  • 


(2.26) 


By  total  differentiation# 


ut  ■  -Ctt*  *  nt  “  'cnx  • 


(2.27) 


where  C  is  the  speed  of  the  wave  front*  From  (2*24 ),  (2*25)  and  (2*26),  it  follows  that 


ut  *  \  *  0  *  \  *  UXA  /B"  •  0 


Comparing  (2*2?)  and  (2.28)#  we  have 


u“  -  c-1l£  ,  c  -  (a‘/b“)1/2  . 


(2.28) 


(2.29) 


Now  we  differentiate  (2*24)  with  respect  to  t  and  (2*25)  with  respect  to  x#  and 
evaluate  the  equations  behind  the  wave  front*  When  we  eliminate  and  make  use  of  the 


expression 


to  obtain 


Hence# 


2  -  *  j  - 

c  u _  -  ufc.  -  c  d(u  )/dx  -  cdu  /dx 

XX  tt  X  t 


-2cd(n[|)d*'1  ♦  <r£)-1  [c'  -  imb'oj  ♦  3c_1  -  o  . 


r  -  /  hx(h2  ♦  h2)‘,/2d*  . 


XX  X 

•0°‘1/2[(f)*0  /  C"5/2**P  f  r“(2A-)~1dx'  ♦  /  r"(2A-)_1dx  ,  (2.20) 


v  '-".vv.  ■r.  ■ 

-V  lA  *  ->  AaV-W 


ifhtrt  «q  is  ths  initial  vslus  of  at  x  -  Xq.  W«  call  x  «  A  a  shoreline  if 

A~(A)  «  0  but  B "(A)  #  0,  and  let 

1 

x  x 

Kx)  *  (j)  /  c"5^2exp  /  r’(2A’)-1dx  • 

XA  XA 

0  0 

Suppose  ag  <  0*  If  1(A)  *  •#  then  \  9  and  the  wave  breaks  before  it  reaches  the 

shoreline*  If  1(A)  *  «•,  then  either  the  wave  breaks  before  it  reaches  the  shoreline  or 

it  breaks  at  the  shoreline*  Next  suppose  ag  >  0*  If  1(A)  *  then  the  wave  breaks  at 
the  shoreline*  Otherwise  if  1(A)  »  •,  we  evaluate  the  lieit  of  nx  given  by  <2*30)  as 
x  ♦  A  and  obtain 

li*  -  (|)  [-<<f)V4  ♦  r“/(2B*)]  ,  (2.31) 

x+A  X 

where  d~  »  A~/B~ «  Hence  the  wave  will  never  break  if  (d  )*  is  finite  at  x  -  A. 

However  for  channels  of  variable  cross  section  the  equilibrium  water  surface  may  converge 

to  a  point  and  this  case  is  also  of  interest*  Assume  again  Sq  >  0,  1(A)  *  •*  Of 
B“(A)  -  d"(A)  -  0,  (d“)’  is  finite  at  x  -  A,  and  h(x,y,z)  *  -z  ♦  g(x,y)  -  0,  we  have 

b2 

r"  -  /  h  <h2  ♦  h2f  1/2d*  -  /  g  dy  , 

L  X  y  *  -bt 

where  y  *  -b^bj  are  the  endpoints  of  the  width  B~(x).  It  follows  from  (2*31)  that 

U‘  \  “  (|)  *  V2]x-t  ' 


and  the  wave  will  never  break 


^  --  V. 


l«  K-dV  KXJATIOH  AMP  THK  INFINITE  HASS  PlUgHA 

Wt  only  sketch  the  derivation  of  the  K-dF  equation*  for  e  channel  of  variable  cross 
section  and  the  details  nay  be  found  in  Zhong  and  Shan  (1982).  we  first  introduce  the  non- 
dinensional  variables 

t  -  0"3/2t*/<H/g)1/2,  (x»y,*>  -  (B_3/2x*/H,y»/H,«*/H)  » 

H|hfp  and  (u,v,w)  are  the  sane  as  before*  The  nethod  used  here  is  a  specialisation  of 
the  procedure  developed  by  Shan  and  Keller  (1973)*  Assune  that  u,v,w,p  and  n  as 
functions  of  t,x,y  and  s  also  depend  explicitly  upon  a  new  variable 

£  -  $S(t,x)  , 

where  8,  a  function  of  t  and  x  only,  will  be  a  phaee  function,  and  that  they  possess 
an  asynptotlc  expansion  of  the  fore 

+(£,t,x,y,s,6)  ~  ♦  0  ♦  $  292  ♦  ...  . 

The  xeroth  approxination  is  assuned  to  be  given  by 

'  <VW  *  °'  *0  *  %  m  0  # 

The  equations  for  the  first  approxination  determine  a  Hamilton -Jacobi  equation  for  S.  Let 

k  -  S  ,  w  -  -S.  . 
x  t 

Then 

0)  *  kG(x) ,  G(x)  -  ±(a(x)/b(x)]  ,  (3.1) 

where  a(x)  is  the  area  of  a  cross  section  Dq,  and  b(x)  is  the  width  of  DQ,  of  water 
at  rest  (Fig.  1).  (3.1)  nay  be  solved  by  the  nethod  of  characteristics  and  the 

corresponding  characteristics  aquations  (Courant  and  Hilbert,  1962)  are 

dt/do  *  M  dx/do  ■  uG(x),  dk/do  -  -k|iG'(x),  dw/d o  •  d S/do  -  0  .  (3.2) 

Where  v  Is  a  proportionality  factor,  we  choose  \i  •  1  so  that  o  -  t.  The  eolutlons  of 
(3.2)  determine  a  one-parameter  family  of  bicharacteristics  called  rays, 

x  *  x(t,o1> 

where  is  constant  along  a  ray.  The  equations  for  the  second  approximation  determine  a 

K-dV  equation  with  variable  coefficients 

Vit  ♦  Vix  ♦  -2ni +  Wic  +  Vu«  ■ 0  *  <3,3) 

Here 


-9- 


■0  ■  2b(x)  ,  (3.4) 

■,  -  2*(x)/G(x)  ,  (3.5) 

■  -  -[G(x)]"1  /  h  (h2  ♦  h2)"1/2<U  -  G~2(x)G* (x)a(x)  ,  (3.6) 

2  x  y  x 

Lo 

■3  -  3k[G(x)]-1b(x)  -  «*,Ui(t,*,r2,0)  -  Uytt.x.y^O)]  ,  (3.7) 

».  »  of1  //  (V#)2dydx  ,  (3.8) 

D0 


Lq  is  the  wetted  boundary  of  Dg#  y  -  y1#y2  are  the  endpoints  of  the  width  of  Dq,  and 

+  is  the  solution  of  the  following  Neuaann  problea 

V2  +  -  k2  in  D0  , 

2 

4  *  w  at  z  *  0  , 

* 

♦  h+  #h-0atL_. 

Ty  y  Tz  *  0 

Since  froa  (3*2) 

d/d a  -  3/3t  ♦  G(x)  3/3x,  dx/do  -  G(x)  , 
along  a  ray#  we  aay  express  ( 3*3 )  in  teras  of  0  and  € 

Vio  +  Vi  +  Vinu  +  View  ■  0  ' 

or  in  teras  of  x  and  £ 

Vix  +  *2ni  +  VtVVi(«‘°  •  (3*9) 

TO  be  definite#  we  choose 

G(x)  -  ♦[*<x)/b(x))1/2  , 

S  *  -t  ♦  /  [G( z) ]  1d*  , 

0 

which  is  a  solution  of  (3.2)  and  It  follows  that 

w  •  1#  k  -  G  \x)  • 

fie  note  that  other  choices  of  S  are  also  possible*  For  rectangular  and  triangular 
channels#  the  coefficients  given  in  (3*4)  to  (3*8)  can  be  explicitly  evaluated  (Zhong  and 
Shen#  1982). 
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Let  d(x)  and  b(x)  be  the  variable  depth  and  width  respectively 

*0  -  2b ( x) ,  -  2b(x>d1/2(x),  ^  -  b'(x)d1/2<x)  ♦  d*  (x)b<x)d"1/2(x)/2  , 

■j  -  3d_1(x)b(x),  »4  -  (1/3)b(x)d(x)  . 

2)  Triangular  channel 

Let  the  two  aides  of  a  cross  section  Dq  be  defined  by  z  -  -  d(x)# 

s  -  -u2(x)y  -  d(x)#  where  nA<x)  -  dCxi/b^x),  i  «  1,2. 

■0  -  2[b1(x)  ♦  b2(x)J  -  2b(x>,  -  G  d,/2(x)b(x)  , 

■j  -  /2  d"1/2(x)[b'(x)d(x)  ♦  d'(x)b(x)/2]/2,  -  5d_1(x)b(x)  , 

■4  -  Id_1(x)/4) [b(x)d2(x)  ♦  <b*(x)  ♦  b2(x)/3]  . 


In  both  cases#  ■1  >  0#  *3  >  0#  *4  >  0  if  d(x)  >  0,  b(x)  >  0*  For  general  cases#  as 

seen  froai  (3.5),  (3*8)  >  0#  m4  >  0  where  we  choose  G(x)  >  0f  but  the  eign  of  *3 

given  by  (3.7)  is  not  obvious.  We  shall  assume  that  *3  *  0  is  the  case. 

Let 


(3.10) 


In  teras  of 


subject  to 


A  -  *3(x)«^1(x)n1  . 

T  and  A#  (3.9)  beccams 

At  ♦  AA ^  +  A^^  -  H(T)A,  T  >  0#  -•  <  J  <  •  , 


(3.11) 


(3.12) 


A(0#S>  -  Aq (  €)#  -•<«<•,  (3.13) 

where 

h(t)  -  -a^xjm^tx)  -  ^ (x)e~2(x)«3(x)  [m4(x)/*3(x)] 1  . 

A  global  existence  theorem  for  (3.12)  and  (3.13)  can  be  easily  obtained  by  extending  the 
existence  results  for  the  K-dV  equation  with  constant  coefficients  due  to  Rato  (1975, 

urns  H(  t)  is  continuous  and  let  H* (-«*#•)  denote  the  Sobolev  space  of  order 


1980).  Assi 


■  of  the  L  -type.  81nee  H9x)A  satisfies  the  conditions  (fl),  (f2)  in  Kato  ( 1975) #  we 
have  the  following  local  existence  result  (Kato,  1980). 

Theorem  1. 


(3.12)  subject  to  (3.13)  with  AQ  S  H8,  s  >  2  has  a  unique  solution 
A  €  C[0,T*  »H*]  c’lOrT'lH*"3] ,  A(0,C)  “*0<5)  , 

for  sone  T*  >  0  and  A(t)  depends  upon  Aq  continuously  in  the  Hs-norm. 

Hereafter  we  shall  denote  an  H*-nor»  by  !•!  and  an  L2-norm  by  1*1.  To  show  that 
there  exists  a  global  solution  in  [0,T]  for  any  T  >  0,  we  need  a  regularity  result  and 
an  a  priori  estinate  for  IA*2,  which  nay  be  obtained  by  naans  of  the  first  three 
conservation  laws  for  the  K-dV  equation  with  constant  coefficients  (Lions,  1969). 

Theoren  2. 

If  A  €  C[0,T>H*]  is  a  solution  to  (3.12)  with  s  >  2  and  if  A(0)  -  AQ  e  H8 
with  s'  >  s,  then  A  6  C(0,T>H8  ]  with  the  sane  T. 

This  theoren  is  a  slnple  extension  of  Kato's  result  (1980)  and  we  omit  the  proof. 


2 

Suppose  for  any  T  >  0  and  t  in  [0,T] ,  a  solution  of  (3.12)  A  €  C[0,T;H  ] 

2 

exists  and  A(0)  *  AQ  €  H  •  Then 

■A(t)I2  <  •  (IA(0)I2)  , 

where  +2  T< • )  depending  upon  T  is  a  nonotone  increasing  function  with  *2  T(0)  *  0. 
Proof i  Assism  first  A,Aq  €  H4.  We  Multiply  (3.12)  by  A  and  integrate  with  respect  to 
£  fron  to  •  to  obtain,  for  t  in  (0,T], 

00 

dlAl2/dt  -  2  /  H(T)A2dC  . 


It  follows  that 


and 


IAI  <  ( IA(0) I ) 


(3.14) 


♦.  (IA(0)I)  -  IA  laxp(MT)  , 

0,1  0 

M  -  SUp  |H(  T)  |  • 

0<T<T 

Haxt  v*  multiply  (3.12)  by  A2  ♦  2A^  and  intagrata  again  to  obtain 

m  m 

d/dT  /  <A3/3  -  A2)d$  -  H(t)  /  (A3  -  2Ag)d(  . 

Than  «•  intagrata  tha  abova  aquation  with  raspact  to  t  fro*  0  to  x  to  obtain 
IA  I2  <  lA.(O)  I2  ♦  (1/3)  IA(0)  I3  +  ( 1/3)  lAl3 

^  *  T.  T. 


Sinca 


1  m  7  T  m  % 

*  2  }  | H ( T* )  1  /  AgdCdT*  +  /  |H(T')|  /  |Ard5dt'  . 

0  ••  ’  0  -• 


IAI33  <  2IA^I1/2IAI5/2  <  lA^I2  ♦  C0 
L 


5/2 

whara  C_  -  2#,. '  #  wa  hava 
0  0#T 


IA?»2  <  (j)  IA?(0)I2  ♦  (|)  IA(0)I5/2  +  C0/3  ♦  «C0  *  3H  }  U^/dT* 


It  follows  that 


,A,1  *  C1.T(,*0,1>  *  ^  "  f  ,A,?dT'  » 

whara  is  a  monotom  incraasing  function  with  C^T(0)  -  0.  By  Gronwallvs 

inequality  , 

IAI,  <  C^2nxp(9*T/4)  -  *t  (»A*1)  .  (3.15) 

Similarly,  «a  multiply  (3.12)  by  A3  ♦  3A^  ♦  6AA^  ♦  { 18/5)A__„,  intagrata  tha  raaulting 
aquation  and  uka  usa  of 

lAlf  <  2IA.IIAI3  , 


and  0.14),  0.15)  to  obtain 


IM*  <  lA(O)  I*  ) 


For  A#A0  «  H2,  we  may  approximate  Aq,  A  by  sequences  of  smoother  functions ,  make  uaa  of 
tha  result  of  local  continuous  dependence ,  and  theorem  4  to  oomplete  the  proof.  By  Lemma 
1,  Theorems  1  and  2  we  have  the  global  existence  theorem. 

Theorem  3.  For  any  T  >  0,  (3.12)  possesses  a  unique  solution 

A  e  C[0,TjH*]  n  c'lOfTiH*"3]  satisfying  A(0,C>  -  AQ  €  H*,  s  >  2.  A  depends  upon  Aq 
continuously  in  H*  norm. 

2 

It  Is  evident  from  the  above  theorem  that  if  we  prescribe  A^  *  a  each  0£  at  x  -  0,  then 

lim  A  -  0  for  any  x  >  0.  If  there  is  a  shelf  created  behind  the  solitary  wave,  it  can 
€+*• 

never  be  extended  to  infinity  in  finite  t. 
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4.  FlgglOW  or  80LIT0M8 

Nt  consider  an  opan  channal  with  two  uniform  sections  SQ  for  x  <  0  and  S1  for 
x  >  connected  by  a  transition  section  for  0  <  x  <  *1#  and  assume  that  the 
transitions  at  x  -  0  and  x  •  x1  are  sufficiently  smooth  so  that  Theory  3  holds,  and 
that  Xj  >  0  is  sufficiently  small*  Let 


-1 

1  -  m(x)B#  m(x)  -  exp  -  /  m  (x§ )m  (x')dx’  . 

O 


In  terms  of  B#  (3*9)  becomes 


B^  ♦  P(x)BBc  ♦  (1/6)Q(s)BU5  -  0  # 


P(x)  -  (x)m(x)m^  (*)#  (1/6)Q(x)  -  {*)  •  (4*3) 

In  81  and  S2t  P  and  Q  are  oonstants  and  we  may  transform  (4*2)  to  the  standard  form 
of  the  K-dV  equation  with  constant  coefficients  except  the  change  of  a  sign  (Muira  et  al.  # 


1974).  Let 


0  -  PBA)*  C  -  Q*/*  #- 


and  (4*2),  in  terms  of  a  and  becomes 

♦  6UU^  ♦  0^  -  0  ,  (4.5) 

Suppose  we  envisage  a  solitary  wave  moving  from  x  *  -•  toward  x  ■  as  a  progressive 
solution  of  (4.2)  in  81#  and  would  like  to  find  the  condition  under  which  fission  of 
soli  tons  will  take  place  after  the  solitary  wave  moves  into  s2. 

A  progressive  wave  solution  of  (4.5)  is  given  by 

0  -  (c/2)MCh2[(e)1/2(C  -  c«>/2]  .  (4.6) 

In  othr  words#  if  we  prescribe  U  at  ;  *  0  in  the  form  of  (4.6)#  0  will  be  a  solution 


of  (4*5)  for  all  C«  Assume 


Prom  (4.1)#  we  have 


B  -  Bq  sech  o£  at  x  -  0 


n1  *  Bgm(0)sech  at  x  *  0  # 
where  m(x)  -  m(0)  for  x  <  0.  Then  (4.4)  implies 
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(4.8) 


0  -  (P< 0  )B0/QC 0 )  J  ««ch2aC  at  C  -  0  , 
where  P(x)  -  P(0),  Q(x)  •  Q(0)  for  x  <  0.  In  comparing  (4*8)  with  (4*6),  Aq  and  a 
mast  satisfy  ths  following  conditions 

1/2 

P(0 )Bq/Q0  -  c/2,  c  ' /2  -  o  , 

so  that  (4*7)  is  ths  initial  condition  for  a  progressiva  solution  of  (4.5) •  Hence  it 
follows  that 

P(0)B0/[2Q(0)«2)  -  1  .  (4.9) 

After  the  solitary  wave  novas  into  S2,  we  neglect  the  effect  of  the  transition 
section.  With  (4.4)  given  by 

0  -  P(x1  )B/Q(x1 ) ,  C  *  Qfx^x/6  ,  (4.10) 

where  P9x)  -  P(x1)r  Q(x)  -  Q(x<|)  for  x  >  x^,  (4.5)  nay  be  solved  by  the  Inverse 

scattering  method,  and  the  corresponding  eigenvalue  problem  is  posed  by  the  equation 

d 2*/dC2  ♦  (U(0,C)  ♦  X)t  -  0,  -•<*<•#  (4.11) 

where  by  (4.4)  and  (4.7), 

0(0, C)  -  (P(Xj )Bq/Q(Xj )  J  sech2a£  •  (4.12) 

From  the  known  results  (Tappert  and  Zabusky  1971;  Johnson  1973i  Muira  at  al.  1974),  U 

consists  of  n  soli tons  for  large  c  if 

P(x1)B0/tQ(*1)«2]  -  n(n  +1),  n  -  1,2,...  (4.13) 

2  2 

Xm/a  ■  -a  ,  x  ■  1,2,...,n  • 
and  the  amplitude  of  the  soli ton  is  given  by 

•  -21^,  m  -  1,2,  *..,n  • 

We  may  eliminate  Bg/o2  from  (4.9)  and  (4.13)  to  obtain 

P(x  )Q(0) 

p(0)Q(x  }  -  n(n  +  1 )/2  n  -  1,2,...  ,  (4.14) 

which  yields  a  relationship  among  P(0),  P(x1),  Q(0),  Q(x^)  for  the  fission  of  n 
soli tons  to  take  place.  For  rectangular  and  triangular  channels,  P(x)  and  Q(x)  can  be 
explicitly  expressed  in  terms  of  channel  width  b(x)  and  depth  d(x),  and  (4.15)  assumes 
rather  simple  forms. 


Let  d(x>,  b(x)  *11  um  constant  values  d^,  bg  in  S0  and  d1#  b1  in  S1»  It  ia 
found  that 

■(x)  -  b"1/2(x)d"1/4<x>,  P(x)  -  (|)d"1/4<x)b"1/2<x),  (j)Q(x)  -  (£)d-1/2  . 

(4.14)  now  becomes 

(d1/d0)"9/4(b1/b0)"1/2  -  n(n  +  1)/2,  n  -  1.2,...  .  <4.15) 

If  b,  -  bQ,  (4.15)  r*doo«s  to  the  criterion  obtained  by  Tappert  and  Zabusky  (1971),  and 
Johnson  (1973). 

(2)  Sy— trie  triangular  channel# 

Let  the  two  sides  of  the  channel  be  defined  by  t  -  ±w(x)y  -  d(x),  where 
ti(x)  -  2d(x)/b(x>  and  let  d(x),  b(x)  assume  constant  values  dg,  b0  and  dv  bt  in 
8q  and  respectively e  Then  (4e14)  reduces  to  (Zhong  and  Shen,  1982) 

[(Y2  ♦  12)/(t^  ♦  12)l(Y0/Y1)"1/2<«1/^)r11/4  -  n(n  ♦  1)/2,  n  -  1,2,...  , 
where  ya  *  b^/d^  A 

In  solving  the  K-dV  equation  (4.5)  with  oonstant  coefficients  by  the  inverse 
scattering  method,  we  tacitly  assumed  that  the  Initial  condition  used  in  (4.11)  was  given 
by  (4.12)  and  the  effect  of  the  transition  section  upon  the  evolution  of  the  solitary  wave 
was  completely  neglected.  However,  if  is  sufficiently  small,  then  by  the  continuity 

of  h  and  its  continuous  dependence  on  Aq  as  obtained  in  Theorem  3,  the  approximation 
made  here  is  justified.  Let  n^x^C)  be  the  solution  of  (3.9)  at  x  -  x1  with  the 
initial  condition  -  ^(0,5)  «  H^,  s  >  2  and  n*(x,C)  and  n^x,*)  are  respectively 
the  solutions  of  (3.9)  for  x  >  x^  with  initial  conditions  nJfx^C)  *  ^(0,0, 
n1  *  that  the  co***1^**1**  ln  (3.3)  are  sufficiently  smooth.  We  state 

the  result  as 
Theorem  4. 

For  any  given  X  >  0  and  c  >  0,  there  exists  5  >  0  such  that 

inf(x,5)  -  <  C  for  x1  <  a,  X  e  [0,X],  and  s  >  2  . 

He  note  that  the  theorem  cannot  tell  what  happens  as  x  ♦  •#  This  concerns  the  global 


stability  of  a  solution  of  the  K-dV  equation,  which  is  an  interesting  problem  for  further 
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On  the  basis  of  some  approximate  equations  derived  for  nonlinear  water 
waves  in  a  channel  of  variable  cross  section ,  several  problems  of  interest , 
which  have  attracted  much  attention  recent ly,  will  be  taken  up  in  this 
report.  First  we  shall  consider  the  breaking  of  an  acceleration  wave  moving 
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channel  of  variable  cross  section,  will  be  resolved.  Finally  we  shall  use  an 
approximation  method  to  study  the  fission  of  solitons  in  a  general  channel  and 
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